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Abstract 

We review two different theoretical approaches to the strong interaction dynamics at the early 
times immediately following heavy ion collisions. One approach is based on small-coupling 
physics of the Color Glass Condensate (CGC). The other approach is based on Anti-de Sitter 
space/Conformal Field Theory (AdS/CFT) correspondence and may be applicable to describ- 
ing large-coupling QCD interactions. We point out that in terms of theoretical tools the two 
approaches are somewhat similar: in CGC one deals with classical gluon fields produced in a 
nuclear shock wave collision, while in AdS/CFT one studies classical gravity in a gravitational 
shock wave collision. We stress, however, that the resulting physics is different: the classical 
gluon fields in CGC lead to a free-streaming medium produced in heavy ion collisions, while 
the classical gravity in the 5-dimensional AdS bulk is likely to lead to ideal hydrodynamics de- 
scription of the produced medium. Also, the valence quarks in colliding nuclei in CGC continue 
along their light cone trajectories after the collision with very little recoil, while we show that in 
AdS the colliding nuclei are likely to lose most of their energy in the collision and stop. 



1. Classical Gluon Fields in CGC 

We begin by considering a high energy heavy ion collision. The basic premise of the CGC 
physics is that for each of the colliding ultrarelativistic nuclei the small-x wave function is char- 
acterized by the large transverse momentum scale called the saturation scale and denoted by Q s . 
At high energy and/or for large nuclei the saturation scale is large, Q s » Aqcd, making the 
strong coupling constant small, a s <k 1, thus allowing for a small-coupling description of the 
small-* nuclear wave functions. Collisions of two nuclei with such wave functions would lead to 
interactions characterized by perturbatively large saturation scales as well, allowing for a small- 
coupling description of the early stages of heavy ion collisions. We refer the reader to HI for a 
review of CGC physics in nuclear collisions. 

If one is interested in the dynamics of the medium produced in the collisions over a not very 
broad rapidity interval (Ay < 1 /a s ), then the relevant collision dynamics is described in the 
framework of the McLerran-Venugopalan (MV) model @]. The MV model states that, due to 
the high parton density in the colliding nuclei, the dominant gluon fields produced in heavy ion 
collisions are classical, and are described by the classical Yang-Mills equations 

D V F" V = P (1) 



with the source current J 11 given by the color charges in the colliding nuclei. This setup is 
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Figure 1: Heavy ion collision in McLerran-Venugopalan model. 

illustrated in Fig. Q] where the two nuclei moving away from each other provide the source 
current, and the classical gluon field is left behind the two nuclei. 

The exact analytical solution of Yang-Mills equations (Q]) providing gluon field generated in 
a heavy ion collision does not exist due to the complexity of the problem. In the diagrammatic 
language to find the classical gluon field one has to resum an infinite set of Feynman diagrams, 
and example of which is shown in Fig. [2] There exist however perturbative solutions |3i|4jj5|], an 
analytic solution for the gluon production cross section in proton-nucleus (pA) collisions [6] and 
a numerical solution of the full nucleus-nucleus (AA) problem [7]. There is also an analytical 
ansatz for the full solution for the gluon production cross section in AA collisions [8|]. 




Figure 2: Diagrammatic representation of the classical gluon field produced in a heavy ion collision in McLerran- 
Venugopalan model. The thick crosses represent nucleons in the nuclei with the top row of the crosses denoting one 
nucleus and the bottom row denoting the other nucleus. The thin cross denotes the point in space-time where we measure 
the gluon field. Incidentally, the same diagram describes classical graviton field produced in a collision of two shock 
waves in AdSs . 

To understand the matter produced by the classical gluon fields in AA collisions let us first 
note that the distribution of this classical matter is rapidity-independent. In the MV model the 
nuclei have a very large transverse extent and are translationally invariant in the transverse di- 
rection. The matter distribution thus does not depend on the transverse coordinate x± and on 
the space-time rapidity rj = (l/2)ln[(x° + x 3 )/(x° - x 3 )], and depends only on the proper time 
r = V(* ) 2 - O 3 ) 2 - (Here x° is time and x 3 is the collision axis: see Fig.[3]forthe explanation of 
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Figure 3: The space-time picture of the ultrarelativistic heavy ion collision in the center-of-mass frame. The collision 
axis is labeled x 3 , the time is x°. 



the coordinates.) The most general energy-momentum tensor of a matter distribution dependent 
only on r can be shown to be of the following form at mid-rapidity (x 3 = 0) |9s] 
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(2) 



in the x°,x l ,x 2 ,x 3 coordinates. In general the transverse pressure p(f) is not equal to the lon- 
gitudinal pressure pi(r). Imposing the energy-momentum conservation condition d^T^ — on 
Eq. <£3J> yields 

£ = -i±£. 0) 

CIT T 

Classical gluon field dynamics is conformal: therefore T% = which implies e - 2p + p^,. 

At early proper times t «: 1 /Q s the classical gluon fields lead to the following scaling of 
energy density e(r) II 1 Qi 1 1 ill 

e(r) ~ ln2 4r- ( 4 ) 

t«\/Q, T ^ s 

With the help of Eq. (01 this leads to the following early-time energy -momentum tensor of the 
produced matter 

e(r) ^ 

e(r) 

e(r) 

-e(r) 



pi" 



T«\/Q, 



(5) 



Note that the longitudinal pressure at early times is negative pi = -e. This result is true for any 
rapidity-independent medium distribution at early time which has a finite total energy. 
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At late proper times t » 1/2? both the analytical perturbative approaches js}] and the full 
numerical simulations [7] lead to the energy density scaling as 



T»l/Q, 

This gives the following energy-momentum tensor 



77" 



t»1/Q, 



(6) 
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(7) 



This late-time energy-momentum tensor has zero longitudinal pressure and corresponds to free- 
streaming of the produced medium. One can see that by arguing that the net energy of the 
produced medium is E ~ ex. Hence E is constant for the energy density from Eq. (O, which can 
be understood as being due to non-interacting particles free-streaming away from the collision 
point. 



2. Classical Gravity in AdSs 



The classical CGC picture of heavy ion collisions is self-consistent, in the sense that it as- 
sumes that a s <k 1 and then justifies the assumption by generating a large momentum scale Q s . 
However, it lacks an essential ingredient needed to describe heavy ion collisions: it does not lead 
to ideal hydrodynamics, which is known to describe RHIC data on particle spectra and elliptic 
flow rather well |Q~2, 13]. In the rapidity-independent case considered above, ideal hydrodynam- 
ics was first considered by Bjorken 01411 . It has 
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(8) 



For the ideal gas equation of state e = 3p it has 



(9) 



Obtaining (ideal) hydrodynamics in a first-principles calculation for heavy ion collisions is the 
important problem of thermalization and/or isotropization of the produced medium. Below we 
discuss this problem assuming that strong-coupling non-perturbative QCD effects are responsible 
for the onset of the hydrodynamic behavior. Since it is unknown how to consistently study QCD 
at strong coupling using analytic methods we will instead use AdS/CFT correspondence, which 
would allow us to study a distant cousin of QCD — the N — 4 super- Yang-Mills (SYM) theory. 

AdS/CFT correspondence conjectures that the dynamics of N — 4 S U{N C ) SYM theory in 
four space-time dimensions is dual to the type IIB superstring theory on AdSsxSs II 1 511 - In the 
limit of large number of colors N c and large 't Hooft coupling A = g 2 N c (with g the gauge cou- 
pling constant) such that N c » A » 1, AdS/CFT correspondence reduced to the gauge-gravity 
duality: N = 4 SU(N C ) SYM theory at N c » A » 1 is dual to (weakly coupled) classical 
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supergravity in AdSs. Hence the gauge dynamics at strong coupling, which includes all-orders 
quantum effects, is equivalent to the classical dynamics of supergravity. Instead of summing 
infinite classes of Feynman diagrams in the gauge theory or using other non-perturbative meth- 
ods, one can simply study classical supergravity in 5 dimensions. For a review of AdS/CFT 
correspondence see fl^ll . 

In a real-life heavy ion collision at RHIC the early-time dynamics is likely dominated by 
the weak-coupling CGC physics. Strongly-coupled dynamics may set in only at later times 
t ~ 1 /Aqcd, though even this time estimate is very crude. As can be shown using the techniques 



of 11711 . matching of perturbative CGC physics onto AdS/CFT dynamics at later times is not 
unique, and does not lead to a single uniquely defined dual geometry in AdSs, allowing instead 
for a variety of possible metrics. Indeed the expectation value of a single local operator (the 
energy momentum tensor) coming from CGC can not uniquely constrain the full quantum state 
of the field theory. To keep our calculations under theoretical control we will consider the whole 
collision of two nuclei in the strong coupling AdS/CFT framework, understanding that this is 
simply a rough approximation of the real heavy ion collision. 

Our goal is to describe the isotropization (and thermalization) of the medium created in heavy 
ion collisions assuming that the medium is strongly coupled and using AdS/CFT correspondence 
to study its dynamics. We want to construct a metric in AdSs which is dual to an ultrarelativistic 
heavy ion collision as pictured in Fig. [3] 



We start with a metric for a single shock wave moving along a light cone [18]: 

ds 2 = ^ l-2dx + dx~ + <r__0O>z 4 dx- 2 + dx\ + cfe 2 j . (10) 

Here x ± = Z is the coordinate describing the 5th dimension such that the boundary of the 

AdS space is at z = 0, and L is the curvature radius of the AdS space. According to holographic 

renormalization 11911 . (T (x~)) is the expectation value of the energy-momentum tensor for a 

single ultrarelativistic nucleus moving along the light-cone in x + -direction in the gauge theory. 



Figure 4: A representation of the metric 1101 as a graviton (wavy line) exchange between the nucleus at the boundary of 
AdS space (the solid line) and the point in the bulk where the metric is measured (denoted by a cross). 

The metric in Eq. (TlOb is an exact solution of Einstein equations in AdSs: R^y + j$ g MV = 
0. It can also be represented perturbatively as a single graviton exchange between the source 
nucleus near the AdS boundary and the location in the bulk where we measure the metric/graviton 
field. This is shown in Fig. H] where the solid line represents the nucleus and the wavy line 
is the graviton propagator. Incidentally a single graviton exchange, while being a first-order 
perturbation of the empty AdS space, is also an exact solution of Einstein equations. This means 
higher order tree-level graviton diagrams are zero (cf. classical gluon field of a single nucleus in 
covariant gauge in the CGC formalism B20I1 ). 
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Figure 5: Diagrammatic representation of the metric in Eq. jilt . Wavy lines are graviton propagators between the 
boundary of the AdS space and the bulk. Graphs A and B correspond to the metrics of the first and the second nucleus 
correspondingly. Diagram C is an example of the higher order graviton exchange corrections. 



Now let us try to find the geometry dual to a collision of two shock waves with the metrics 
like that in Eq. d). We will follow jH. Defining h(x~) = ^ {Ti __(*-)> and t 2 {x + ) = 



2jf {T2++{x + )) we write the metric resulting from such a collision as 
L 2 [ 

ds 2 = —I - 2dx + dx~ + dx 2 ± + dz 2 + ti(x~)z 4 dx~ 2 + t 2 (x + )z 4 dx +2 



z 2 



+ higher order graviton exchanges j (11) 

The metric of Eq. (fTTT l is illustrated in Fig. [5] The first two terms in Fig. [5] (diagrams A and 
B) correspond to one-graviton exchanges which constitute the individual metrics of each of the 
nuclei, as shown in Eq. (TTOb . We need to calculate the next order correction to these terms, which 
is shown in the diagram C in Fig. [5] 

Fig.[5]illustrates that construction of dual geometry to a shock wave collision in AdSs consists 
of summing up all tree-level graviton exchange diagrams. It is similar diagrammatically to the 
classical gluon field formed by heavy ion collisions in CGC The diagram for the gluon 

field shown in Fig.|2]can also be interpreted as the diagram for gravitons in AdSs bulk leading to 
the metric produced in a collision of two gravitational shock waves. 



Without going into details of the calculation which can be found in 112 111 one can calculate 
the graph in Fig. [St as follows. Let us take delta-function shapes for the shock waves: t\(x~) = 
/ii S(x~) and t2(x + ) = fi2 S{x + ). The contribution of diagram C in Fig. |5]to the energy density of 
the produced medium e should therefore be proportional to yui fi2- As jUi and jj.2 have dimensions 
of mass cubed each (M 3 ), and dimension of e is M 4 , we need to multiply [i\ ji2 by the square of 
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Figure 6: T ++ component of the proton's energy-momentum tensor after the collision as a function of light cone time x 
(arbitrary units). 



proper time, r , as this is the only dimensionful variable left. One obtains [22, 21] 



e(r) 



(12) 



A more detailed calculation fixes the prefactor in Eq. ( TT2b 1 22 , 2j]] and also specifies the region of 
validity of this result. Energy density is rapidity-independent, because of the cancellation of ra- 
pidity factors due to graviton exchanges in Fig. [3jZT. The energy-momentum tensor corresponding 
to the energy density from Eq. (fT2l) is 
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(13) 



and also has large negative longitudinal pressure, similar to the CGC one (0. There is a problem 
with energy density scaling as shown in Eq. (fT2l (or, equivalently, with the energy-momentum 
tensor in Eq. (fT3l): as was shown in I18I1 . in this case there is a frame characterized by time 
4-direction f in which the energy density is negative f f < 0. Hence the scaling of energy 
density in Eq. ( fT2l with proper time leads to negativity of energy density in other (boosted) 
frames. At this point it is not clear whether this result presents a problem, as there may be 
nothing wrong with energy density becoming negative for a short period of time lf22ll . 

To better understand dynamics of the shock wave collisions let us follow one of the shock 
waves after the interaction. First we "smear" the delta-function profile of that shock wave: 



tl ( x -) = ^-6( X -)6(ai -x~). 

A 1 ' 3 



(14) 



Here p\ oc p + A 2 A 1 ^ and ci\ oc R A oc where the nucleus of radius R has A nucleons in 
it with N}. valence gluons each. p + is the light cone momentum of each nucleon and Aj is the 
typical transverse momentum scale. The "++" component of the energy-momentum tensor of a 
shock wave after the collision at x~ — a\j2 is [21] 



(T ++ (x + ,x- =ai/2)> = ^^\\-2p 2 x +2 ai \. 



(15) 



The first term on the right of Eq. dl~5b is due to the original shock wave while the second term 
describes energy loss due to graviton emission. Eq. dl~5b shows that {T ++ ) of a nucleus becomes 
zero at light-cone times (as p2 °c p ~ A 2 A 1 ^ as p\ in the center-of-mass frame) 

(16, 



^£Ha7 AA 1 / 3 ' 

Zero (T ++ ) would mean stopping of the shock wave and the corresponding nucleus. The re- 
sult can be better understood by doing all-order resummation of graviton exchanges with one 
shock wave, which is needed for modeling proton-nucleus collisions [23]. The full result for the 
proton's "++" component of the energy-momentum tensor is 

Nj. m 1 

< T > ~ T"T / = =. ^ 0<x-<a 1 . (17) 

Eq. ( fTTI i is illustrated in Fig. [6] in which one can see that the proton loses all of its light cone 
momentum over a rather short time. 

We thus conclude that the collision of two nuclei at strong coupling leads to a necessary stop- 
ping of the two nuclei shortly after the collision. If the nuclei stop completely in the collision, the 
strong interactions between them are almost certain to thermalize the system, probably leading 
to Landau hydrodynamics [24]. It is possible that the mid-rapidity region of such collision may 



be well-described by Bjorken hydrodynamics [14], but this still remains to be shown 
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